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Abstract 


The  dynamic  response  of  magnetically  levitated  (maglev)  ground 
transportation  systems  has  important  consequences  for  safety  and  ride  quality, 
guideway  design,  and  system  costs.  Ride  quality  is  determined  by  vehicle 
response  and  by  environmental  factors  such  as  humidity  and  noise.  The  dynamic 
response  of  the  vehicles  is  the  key  element  in  determining  ride  quality,  and 
vehicle  stability  is  an  important  safety-related  element.  To  design  a  proper 
guideway  that  provides  acceptable  ride  quality  in  the  stable  region,  vehicle 
dynamics  must  be  understood.  Furthermore,  the  trade-off  between  guideway 
smoothness  and  the  levitation  and  control  systems  must  be  considered  if  maglev 
systems  are  to  be  economically  feasible.  The  link  between  the  guideway  and  the 
other  maglev  components  is  vehicle  dynamics.  For  a  commercial  maglev  system, 
vehicle  dynamics  must  be  analyzed  and  tested  in  detail.  In  this  study,  the  role  of 
dynamics  and  controls  in  maglev  vehicle/guideway  interactions  is  discussed,  and 
the  literature  on  modeling  the  dynamic  interactions  of  vehicle/guideway  and 
suspension  controls  for  ground  vehicles  is  reviewed.  Particular  emphasis  is 
placed  on  modeling  vehicle/guideway  interactions  and  response  characteristics  of 
maglev  systems  for  a  multicar,  multiload  vehicle  traveling  on  a  single-  or  double¬ 
span  flexible  guideway,  including  coupling  effects  of  vehicle/guideway, 
comparison  of  concentrated  and  distributed  loads,  and  ride  comfort.  Different 
control-law  designs  are  introduced  into  vehicle  suspensions  when  a  simple  two- 
degree-of-freedom  vehicle  model  is  applied.  Active  and  semiactive  control  designs 
for  primary  and  secondary  suspensions  do  improve  the  response  of  vehicle  and 
provide  acceptable  ride  comfort.  Finally,  future  research  associated  with 
dynamics  and  controls  of  vehicle/guideway  systems  is  identified. 


1  Introduction 


A  high-speed  ground  transportation  system,  based  on  magnetically  levitated 
(maglev)  vehicles  propelled  by  a  linear  electric  motor,  has  been  proposed  to  meet 
future  intercity  transportation  requirements.  One  possible  and  attractive 
approach  is  in  replacing  air  travel  for  selected  intercity  trips  of  100  to  600  miles. 
The  maglev  system  will  offer  the  advantages  of  lower  noise  and  emissions  and 
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better  ride  quality,  as  well  as  potential  energy  savings  and  economic  benefits 
(Bohn  and  Steinmetz  1985;  Chen  et  al.  1992;  Coffey  et  al.  1991;  Johnson  et  al.  1989; 
Katz  et  al,  1974;  Zicha  1986). 

While  some  design  concepts  have  been  developed  nearly  to  commercial 
application,  the  attractiveness  of  maglev  systems  is  expected  to  be  enhanced  even 
further  over  the  next  several  years  by  new  or  improved  concepts,  improved  design 
and  construction  methods,  and  new  material  (including  high-temperature 
superconductors,  high-energy  permanent  magnets,  and  advanced  material  for 
guideways).  It  is  therefore  reasonable  to  expect  that  maglev  systems  may  indeed 
be  a  key  transportation  mode  in  the  21st  century  (Chen  et  al.  1992). 

For  several  decades,  research  and  development  have  been  performed  in  the 
areas  of  magnetic  levitation,  response  of  maglev  vehicles  to  rough  guideways, 
interaction  of  variously  suspended  vehicles  with  flexible  guideways,  and 
optimization  of  vehicle  suspensions.  The  results  of  these  efforts  are  useful  in 
providing  appropriate  criteria  for  the  design  of  maglev  systems  (Bohn  and 
Steinmetz  1985;  Chiu  et  al.  1971;  Iguchi  and  Hara  1985;  Katz  et  al.  1974;  Sinha 
1987). 

The  dynamic  response  of  magnetically  levitated  vehicles  is  important 
because  of  safety,  ride  quality,  guideway  design,  and  system  cost.  More  emphasis 
should  be  placed  on  guideway  design,  because  the  cost  of  the  guideway  structure 
is  expected  to  be  60-80%  of  the  overall  initial  capital  investment  cost  (Uher  1989; 
Zicha  1986).  Thus,  guideway  design  is  a  critical  area  of  potential  capital  savings. 
More-flexible  guideways  are  less  expensive,  but  cause  complex  vehicle/guideway 
interactions  and  affect  ride  quality.  An  optimized  guideway  design  will  be 
important  for  a  high-speed  maglev  system  that  offers  good  ride  quality.  As 
maglev  vehicle  speeds  increase  to  200-300  mi/hr,  or  as  guideways  become  lighter 
and  more  flexible  to  reduce  costs,  the  dynamic  interactions  between  vehicle  and 
guideway  become  an  important  problem  and  will  play  a  dominant  role  in 
establishing  vehicle  suspension  requirements  and  specifications  for  guideway 
stiffness,  weight,  and  span  length  (Cai  et  al.  1992a,  1992c;  Chiu  et  al.  1971;  Vu- 
Quoc  and  Olsson  1989;  Zicha  1986). 

Light  guideways,  especially  those  made  of  steel,  may  be  susceptible  to 
dynamic  instability  and  unacceptable  vibration,  and  thus  dynamic  evaluation 
must  be  included  in  the  structural  analysis.  Different  dynamic  responses  of 
coupled  vehicle/guideway  systems  may  be  observed,  including  periodic  oscillation, 
random  vibration,  dynamic  instability,  chaotic  motion,  parametric  resonance, 
combination  resonance,  and  transient  response  (Chen  et  al.  1992). 
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To  design  a  proper  guidew~y  that  provides  acceptable  ride  quality,  the 
dynamic  interaction  of  vehicles  and  guideways  must  be  understood. 
Furthermore,  the  trade-off  between  guideway  smoothness  and  design  of  the 
levitation  and  control  systems  must  be  considered  if  the  maglev  system  is  to  be 
economically  feasible.  The  coupled  vehicle/guideway  dynamics  are  the  link 
between  the  guideway  and  the  other  maglev  components.  Thus,  reliable 
analytical  and  simulation  techniques  are  needed  in  the  design  of  vehicle/ 
guideway  systems  (Cai  et  al.  1992a,  1992c;  Chen  et  al.  1992;  Richardson  and 
Wormley  1974).  Furthermore,  the  coupled  vehicle/guideway  dynamic  model  with 
multiple  cars  and  multiple  loads  must  be  developed  to  meet  the  design 
requirements  of  maglev  systems.  This  analytical  model  should  also  be  easily 
incorporated  into  the  computer  code  for  dynamic  simulation  of  maglev  systems 
(Cai  et  al.  1992a,  1992c). 

Magnetically  suspended  systems  are  intrinsically  underdamped.  Achieving 
safe,  stable  operation  and  acceptable  ride  comfort  requires  some  form  of  vehicle 
motion  control.  Moreover,  vehicle  tolerance  to  guideway  flexibility  and 
roughness,  as  well  as  to  transient  perturbing  forces  such  as  wind  gusts  and 
guideway  misalignments,  will  be  influenced  by  air-gap  size  and  suspension 
control  characteristics,  including  response  time  and  dynamic  range.  To  the 
extent  that  tolerances  can  be  increased  through  suitable  suspension  control 
systems,  guideway  cost  can  be  reduced  and  the  system  made  more  robust  (Cai  et 
al.  1992;  Faye  et  al.  1989;  Kortum  et  al.  1988;  Sinha  1987). 

For  safety,  maglev  systems  should  be  stable.  Thus,  stability  characteristics 
must  be  studied  because  instabilities  in  maglev  system  models  have  been  observed 
at  Argonne  National  Laboratory  and  other  organizations  (Cai  et  al.  1992b;  Chu 
and  Moon  1983;  Moon  1974,  1975).  With  a  better  understanding  of  vehicle  stability 
characteristics,  a  better  control  law  can  be  adopted  to  ensure  a  high  level  of  ride 
comfort  and  safety.  Vehicle  suspension  control  designs  are  therefore  necessary  to 
meet  stability  requirements  of  maglev  systems. 

Although  the  technical  literature  contains  a  substantial  number  of 
publications  dealing  with  controls  of  conventional  ground  vehicles  (Bernard  et  al. 
1987;  Chalasani  1987;  Dukkipati  et  al.  1992;  Elmadany  1990;  Karnopp  and 
Margolis  1984),  very  little  work  exists  on  suspension  control  designs  of  maglev 
systems  (Faye  et  al.  1989;  Gottzein  et  al.  1974;  Katz  et  al.  1974;  Kortum  et  al.  1988; 
Kortum  and  Utzt  1984;  Sinha  1987). 

Therefore,  this  study  is  focused  on  the  dynamics  and  control  of  maglev 
vehicles/guideways.  We  first  discuss  the  problems  associated  with  modeling 
vehicle/guideway  interactions  and  then  explain  the  response  characteristics  of 
maglev  systems  for  a  multicar,  multiload  vehicle  traveling  on  a  single-  or  double- 
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span  flexible  guideway,  with  an  emphasis  on  coupling  effects  of  vehicle/guidewav, 
comparison  of  concentrated  and  distributed  loads,  and  ride  comfort  Second, 
different  control-law  designs  are  introduced  into  vehicle  suspensions  when  a 
simple  two  degree-of- freedom  vehicle  model  is  applied  Active  and  semi  active 
control  designs  for  primary  and  secondary  suspensions  indeed  improve  the 
response  of  vehicle  and  provide  acceptable  ride  comfort  Finally,  future  research 
associated  with  dynamics  and  controls  of  vehicle/guideway  systems  are  identified 


2  Dynamic  Interactions  of  Maglev  Vehicle/Guideway 

Systems _ 

2.1  Background  for  Dynamic  Interactions 

To  simplify  the  vehicle  model,  only  vertical  motions  of  the  vehicle  are 
considered,  based  on  the  assumption  that  vertical  motion  is  dominant  and  that 
other  motions  can  be  ignored  when  vertical  motion  is  evaluated.  This  is 
applicable  in  a  system  in  which  passenger-compartment  vertical  accelerations 
are  limited  to  less  than  0.05  g  and  in  which  vehicle  unsprung  mass  (i.e.,  mass 
associated  with  the  primary  suspension)  inertia  forces  are  low  compared  to 
vehicle  weight.  Thus,  the  influence  of  vehicle  heave  acceleration,  which  is  of 
particular  interest  because  it  is  used  as  a  measure  of  passenger  comfort,  can  be 
determined  in  this  simplified  model  ^Cai  et  al  1992a,  1992c;  Richardson  and 
Wormley  1974). 

In  general,  at  least  a  two-suspension  vehicle  model  is  necessary  to  model 
primary  and  secondary  suspensions  of  maglev  vehicles  (Cai  et  al  1992a,  1992c. 
Richardson  and  Wormley  1974;  Vu-Quoc  and  Olsson  1989 >.  Lumped  masses  and 
passive  parameters,  such  as  linear  springs  and  dashpots,  are  used  to  represent 
these  suspensions.  For  example,  in  most  cases,  the  secondary  suspension  w*as 
described  as  consisting  of  mass,  spring,  and  dashpot,  while  the  primary 
suspension  consists  only  of  spring  or  mass  and  spring  without  dashpot 
(Richardson  and  Wormley  1974).  The  model  developed  in  our  previous  work 
included  masses,  springs,  and  dampings  in  both  primary  and  secondary 
suspensions  (Cai  et  al.  1992a,  1992c).  The  spring  in  the  primary  suspension 
represents  the  magnetic  gap  stiffness,  while  the  damping  represents  a  passive 
damping  control  for  magnets. 

When  vehicle  acceleration  forces  are  much  low’er  than  the  constant  force  due 
to  vehicle  weight,  a  constant  moving  force  can  be  used  to  represent  the  simplest 
vehicle  model.  Both  concentrated  and  distributed  moving  forces  were  used  as 
suspension  forces  in  previous  studies  (Richardson  and  Wormley  1974).  Cai  et  al. 
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(1992a,  1992c)  studied  different  cases  for  various  combinations  of  constant  and 
pulsating  forces  moving  along  the  guideway  and  compared  the  moving-force  and 
quarter-car  models. 

Only  the  single-car  model  can  be  found  in  published  literature  in  analyzing 
vehicle/guideway  interactions.  In  practice,  a  maglev  vehicle  may  include  two  or 
more  cars.  Moreover,  couplings  between  car  bodies  will  certainly  affect  vehicle? 
guideway  dynamics  because  of  constraint  forces.  Therefore,  a  model  that 
represents  a  multiple-car  vehicle  should  be  developed  to  study  dynamic  analysis  of 
vehicle/guideway  systems. 

For  a  flexible  guideway,  elastic  deformation  must  be  considered.  The 
guideway  vertical  motion  is  excited  by  the  full  vehicle  weight,  while  lateral  and 
longitudinal  motions  are  excited  by  only  a  fraction  of  vehicle  weight;  therefore, 
attention  is  focused  on  vertical  guideway  deflection  when  analyzing  vehicle/ 
guideway  interactions.  Beam  theory'  has  been  verified  as  a  good  approximation 
for  guideway  dynamics  when  span  width  to  length  ratio  is  less  than  0.7 
(Richardson  and  Wormley  1974).  Therefore,  the  classical  Bernoulli-Euler  beam 
equation  is  used  to  model  guideway  characteristics  in  virtually  all  recent  analyses 
of  vehicle/guideway  interactions. 

Guideway  surface  irregularities  are  very  important  to  the  ride  quality  of 
maglev  systems.  These  irregularities  may  be  caused  by  imperfections  in 
manufacturing  and  assembling  of  the  structural  components,  as  well  as  by 
thermal  effects  or  surface  wear.  In  practice,  guideway  surface  irregularities  are 
an  important  input  to  maglev  vehicles  and  affect  dynamic  interactions  (Cai  et  al. 
1992a,  1992c).  Guideway  irregularities  can  be  measured  and  statistically  studied 
to  determine  quantitative  relationships  between  tolerances  and  resulting 
guideway  smoothness.  The  guideway  profile  is  the  sum  of  a  static  profile  and  a 
dynamic  profile.  Static  irregularities  depend  on  construction  practice,  settling, 
dead-weight  loads,  and  environmental  conditions  (Snyder  and  Wormley  1977). 
The  dynamic  profile  of  a  guideway  is  composed  of  periodic  and  random  motions  or 
disturbances  due  to  guideway  deflections  from  the  moving  vehicle  iChen  et  al. 
1992). 

Analytical  methods  to  predict  dynamic  vehicle/guideway  interaction  can  be 
divided  into  three  groups:  lumped  mass,  direct  numerical,  and  modal.  The 
lumped  mass  method  is  simple  and  can  be  used  easily  to  account  for  nonuniform 
properties,  while  the  direct  numerical  method  is  accurate  but  requires  more 
computer  time.  The  modal  analysis  method  is  an  efficient  compromise  between 
the  other  two  methods  (Chen  et  al.  1992;  Fryba  1972;  Olsson  1985;  Richardson  and 
Wormley  1974). 
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In  the  modal  analysis  technique,  the  Bemoulli-Euler  equation  is  used  as  the 
basis  for  the  modal  solution  technique  of  distributed  guideway  dynamics  And  the 
space-  and  time-varying  guideway  motion  is  represented  as  an  infinite 
summation  of  the  natural  mode  solutions,  which  is  formulated  as  the  infinite 
sum  of  the  products  of  mode  shapes  and  time-varying  modal  amplitudes  Time- 
varying  modal  amplitudes  depend  on  the  forcing  functions  (i.e.,  interaction  forces 
between  vehicle  and  guideway)  and  initial  conditions. 

In  practice,  a  finite  number  of  modes  are  used  to  represent  guideway  motion. 
The  number  of  modes  required  for  a  given  level  of  accuracy  depends  on  the 
frequency  content  of  the  guideway  forcing  function,  the  traverse  speed  of  this 
function,  and  the  beam  properties. 

Mode  shapes  are  determined  from  the  natural  unforced  vibration  of  the  span 
and  are  affected  by  support  boundary  conditions  at  the  ends  of  the  beams,  and  at 
intermediate  supports  in  the  case  of  multiple-span  beams.  Boundary  conditions 
are  defined  by  the  characteristics  of  the  supports  and  the  coupling  between 
successive  beam  spans.  The  simplest  case  occurs  for  beams  of  only  one  span 
length,  which  are  simply  supported  on  rigid  supports.  If  the  beam  extends  over 
more  than  one  span,  then  at  interior  simple  supports  the  slope  and  bending 
moment  must  be  continuous  across  the  support.  Single-,  multiple-,  and 
continuous-span  guideway  models  for  vehicle/guideway  interaction  in  high-speed 
maglev  systems  can  be  found  in  recent  literature  (Chiu  et  al.  1974;  Smith  et  al. 
1975a,  1975b).  In  most  of  these  studies,  however,  vehicles  were  described  as 
constant  concentrated  or  distributed  force  traveling  along  the  beams.  Cai  et  al 
(1992a,  1992c)  presented  a  detailed  analysis  for  a  two-degree-of-freedom  vehicle 
traveling  along  both  single-  and  double-span  guideways. 

For  maglev  vehicles  restricted  to  vertical  accelerations  of  less  than  0.05  g,  the 
inertia  force  is  much  smaller  than  the  static  load,  normally  about  5 %  or  less. 
Also,  coupling  between  vehicles  and  guideways  will  be  small,  i.e.,  if  the  dynamic 
suspension  forces  acting  on  the  guideway  are  low  compared  with  the  static  force 
due  to  vehicle  weight,  dynamic  coupling  will  be  low  (Richardson  and  Wormley 
1974).  In  this  case,  the  guideway  deflection  profile  is  computed  by  assuming  that 
the  suspension  forces  are  constant  at  their  static  values  and  move  along  the 
guideway  at  vehicle  speed.  The  deflection  is  then  used  as  known  displacement 
input  into  the  suspension,  and  the  vehicle  dynamic  motions  are  determined  by 
standard  transfer  function  analysis. 

When  the  unsprung  mass  is  greater  than  25%  of  the  vehicle  mass,  such  as  in 
an  electromagnetic  system  (EMS),  or  when  vertical  vehicle  accelerations  can  be 
greater  than  0. 1  g,  guideway  deflection  may  be  significantly  affected  by  dynamic 
suspension  forces,  and  fully  coupled  analysis  of  vehicle/guidewa>  interaction  is 


7 


needed.  In  the  EMS  system,  the  large  accelerations  of  the  primary  suspension 
system  mass  due  to  guideway  roughness  cause  significant  excursions  in  magnet 
reaction  force,  and  it  appears  that  the  vehicle/guideway  equations  should  not  be 
decoupled  (Katz  et  al.  1974).  Cai  et  al.  (1992a,  1992c)  extensively  investigated 
coupled  effects  of  vehicle/guideway  interactions  in  a  wide  range  of  vehicle  speeds 
with  various  vehicle  and  guideway  parameters  for  maglev  systems  and  provided 
appropriate  criteria  for  decoupling  at  critical  vehicle  speeds  or  crossing 
frequencies. 

To  evaluate  a  wide  range  of  vehicle  and  guideway  designs  for  an  equally  wide 
range  of  operating  conditions,  it  is  necessary  to  develop  dynamic  models  to 
describe  dynamic  response  of  vehicle/guideway  interactions.  Various  computer 
codes  have  been  developed  to  provide  the  necessary  dynamic  simulations  (Wang  et 
al.  1991).  However,  most  existing  computer  codes  for  maglev  systems  cannot 
simulate  dynamic  interaction  of  vehicle/guideway  because  these  codes  can  only 
model  the  guideway  as  a  rigid  body.  We  recently  incorporated  dynamic 
interaction  of  vehicle/guideway  for  multicar,  multiload  maglev  systems  into  a 
new  computer  program  based  on  the  ANLMAGLEV  program,  which  we 
developed  in  1991  at  Argonne  National  Laboratory.  ANLMAGLEV  can  simulate 
the  dynamic  behavior  of  a  wide  range  of  vehicle  and  guideway  designs  over  a 
broad  range  of  operating  conditions  and  provides  broad  applicability  to  maglev 
system  analysis  (Coffey  et  al.  1991;  Wang  et  al.  1991). 


2.2  The  Vehicle  Model 

A  multicar,  multiload  vehicle  traveling  along  a  flexible  guideway  at  a  velocity 
v,  as  shown  in  Fig.  1,  is  considered  in  our  mathematical  model  for  dynamic 
analysis  of  vehicle/guideway  interactions.  The  car  body  is  rigid  and  has  a 
uniform  mass.  The  center  of  mass  is  consistent  with  that  of  moment  of  inertia. 
Each  car  is  supported  by  certain  numbers  of  magnets  (or  bogies)  with  linear 
springs  and  dampings  (see  Fig.  2),  which  form  the  primary  and  secondary 
suspensions  of  the  vehicle.  If  there  is  only  one  magnet  (i.e.,  the  unsprung  mass) 
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Fig.  1.  Model  of  multicar,  multiload  maglev  vehicle  traveling  along  a  guideway 
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Fig.  2.  Model  of  single  car  supported  with  multiple  magnets  and  traveling  along 
a  guideway 


attached  to  the  vehicle,  there  is  a  single  concentrated  load  and  only  one¬ 
dimensional  motion  (i.e.,  heave  motion)  of  the  vehicle.  If  there  are  multiple 
magnets  on  the  vehicle,  the  loads  are  considered  multiple  or  distributed  and  the 
vehicle  is  capable  of  both  heave  and  pitch  motions.  But  in  this  study,  only  vertical 
vehicle  motion  is  considered  because  it  is  dominant  in  the  dynamic  analysis  of 
vehicle/guideway  interactions. 

The  equations  of  motion  for  the  vehicle  are  then 
N  N 

msySi  +csX(ysi  -  ypij)+ (^si  -yPij) 

j=i  j=i 

+cv{[ySi  -yS(i-i)]+[ySi  -yS(i+i>]} 

(i) 

+  kv{[ysi  -yS(i-i)]+[ySi  -yS(i+i)]}  =  -msg 


(i  =  2, ...,  M-  1;  j  =  1, ....  N) 
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N  N 

msysi  +csX(ysi  -  yPij)+ ks  X(ysi  -yPij) 

j=i  j=i 


+  Cv(ysl-ys2)+kv(ysl-ys2)  =  -msg 


(i  —  1;  j=  h N)  (2) 

N  ^  N 

msySM  +csX(ySM  -yPM.;]+ks  J(ysM  -yPMj) 

j=i  j=i 

+  cv[ysM-ys(M-i)]+kv[ysM -ys(M_1}]  =  -msg  (3) 

(i  =  M;  j  =  1, ....  N) 

and 

mPyPjj +  cP(yPij + y^ij) + kP(yPij + ygij) 

~  rs(ysi  -  yPij)-  ks(ysi  -  ypij) =  ~mP&  ^ 

(i  —  1*  •••»  M;  j  =  I, N) 

where  lumped  masses  mp  and  ms,  linear  springs  kp  and  ks,  and  dampings  cp  and 
cs  represent  primary  and  secondary  suspensions;  the  displacement  of  two 
suspensions  are  yp  and  ys;  subscripts  i  represent  i-th  car  body  and  j  represents 
j-th  magnet  on  the  i-th  car;  M  is  number  of  cars;  N  is  number  of  magnets  on  each 
car;  and  kv  and  cv  are  intercar  stiffness  and  damping,  representing  constraints 
between  adjacent  cars.  For  a  magnetic  primary  suspension,  kp  and  cp  represent 
magnetic  gap  stiffness  and  passive  damping.  ygij  is  guideway  displacement 
input  at  the  i-th  car  and  the  j-th  magnet. 

Uncoupled  natural  frequencies  and  modal  damping  ratios  are  defined  as 
follows: 
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r  =  —IB— 
p  2mpwp  ’ 

r  _  Ncs 
S  2mscos ' 

And  several  nondimensional  parameters  are  introduced: 


Using  Eqs.  5  and  6,  we  can  rewrite  Eqs.  1,  2,  3,  and  4  as 

N  N 

Ysi  "*■  ^Cs^sYsi  ~~  2CS  f  ^  Ypij  +  ^sysi  ~  ^  ^  ^  Ypij 

j=i  j=i 

+  2^scosac[2ysi  ~ys(i-i)  ~y(i+l)] 

+  o)^ak[2ysi  -ys(i_i)  -ys(i+i)]  =  ~g 

(i  =  2, M  -  1;  j  =  1, N) 

N  N 

Ysl  ■*"  ^Cs^sysl  “ 2^scos  / N  ]Typij  +  cogysj  —  0)§  /  N  ]^ypij 

j=l  j=l 

+  2;scosac  (y  s  i  -  y  s2 ) + of  <*k  (y  s  i  -  y  S2 )  =  -e 


(5) 


(6) 


(7) 


(i=  1;  j=l,...,N) 


(8) 
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N 

ysM  +  ^s^s^sM  “  2Cs®s  !  N  ^  y pMi 

j=l 


N 

+«?ySM  -<«>?  /n  XypMj 

j=i 


+  2^sO)sac[ysj^  ys(M-l).  +  c°sak[ysM  ys(M-l)]-  s 


(9) 


and 


(i  =  M;  j  =  1, ....  N) 


y  pij  +  2Cp«P  ( i + Pc  )y  pij + ( 1 + Pk  )y  Pij 

-  2CpO)p(3cysi  -  C0ppkysi  =-g-  2CpC0pygij  -  co^y 

(i  =  1,...,M;  j=  1,...,N).  (10) 

The  system  of  Eqs.  7-10  can  be  represented  in  matrix  form  by  an  M  +  M  x  N  size 
set: 


My  +  Cy  +  Ky  =  Q 


(ID 


where 


Yi  = 


ysi 

yS2 

LysM. 


y2  = 


Mxl 


yPn 

yPi2 

yPiN 

yPMi 

yPM2 

yPMN 


(M-N)xl 


(12) 


12 


Q  = 


[Oil 

[%1 


-el 


L  SjMxl 


~S  -  2^pcopygll  -  ciipYgn 

-g-2;po)pygi2-<flpygi2 


-g-2CpwpygiN-o)2ygiN 


-g  ~  2Cpo)py  gM  i  -  “pygM  l 
~g  —  2CP“PygM2  ~  “pygM2 


y 


2 

-g  -  2CpO)pygMN  ~  wpygMN 


J(MN)xl 


M  = 


C  = 


"1 

1 

o 

C11 

C12 

■C21 

C22 

0 


(M+M-N)x(M+M-N) 


(M+MN)x(M+M  N) 


Cn  = 


2Csws(l  +  ac)  2^g(Oscxc  0  0 

-2CsQ)sac  2Csa)s(l+2ac)  -2Cs«sac  0 

0  -2Cs«sac  2Csa)s(l+2ac)  -2Csa>sac 


c12  ~ 


-2Cs<os/N 

0 


-2CS0)S  /  N 
0 


0 


0 


-2£sa)s/N  -2CS0)S/N 


(12) 

(Contd.) 


IMxM 


Mx(MN) 
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C21  ~ 


-acpcoppc  o 

*  * 

•  • 

-2CpO)ppc  0 

0  -2CpCOppc 

•  • 

•  • 

•  • 

0  -2CpCOpPc 


i{M-N)xM 


C22  ~ 


2^p0)p(l+pc) 

0 

0 


0  0 

2CpC0p(l  +  pc)  0 

0  2Cp(0p(l+pc) 


(MN)x(MN) 


K  = 


*11 

*12 

*21 

*22 

*11  = 


o)f(l+ak) 


'S - -  -“sak 

-®s«k  ©2(l+2ak) 


s 

0 


"®s«k 


0 

"“sak 


0 

0 


©2(l+2ak)  -©2ak 


MxM 


*12  = 


-©s2/n 

0 


-©g/N  0  —  0 

0  -©2/N  -  -co2 /N 

a  S 


Mx(MN) 


*21  = 


-2©pPk 


-2©pPk 

0 

0 


0 

0 

-2©pPk 


‘2©?Pk 


(MN)xM 


(12) 

(Contd.) 
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©J(l+pk)  0  0 

K  _  0  ©Ju+fc)  0 

0  0  ©2(1+  pk)  ... 

:  :  :  :  J(M  N)x(M  N)  (12) 

(Contd.) 

If  M  =  1  and  N  =  1,  the  system  of  Eqs.  11  and  12  will  represent  the  vehicle  to  be 
a  one-dimensional  model  with  two  degrees  of  freedom  (Cai  et  al.  1992a). 


2.3  The  Guideway  Model 


For  typical  guideway  systems,  span-length-to-width  ratios  are  large  enough 
so  that  individual  spans  may  be  considered  as  beams  rather  than  as  plates.  Thus, 
a  Bernoulli-Euler  beam  model  can  be  applied  to  a  freely  supported,  homogeneous, 
isotropic,  and  uniform-cross-section  guideway. 

The  equations  of  motion  for  guideway  spans  where  a  multicar,  multiload 
vehicle  traveling  along  may  be  derived  as 


■d  Yk 


+  m^-^  =  Fk(x,t), 
dt 


where  x  is  the  axial  coordinate  of  the  beams,  t  is  time,  El  is  the  bending  rigidity  of 
the  beams,  C  is  the  viscous  damping  coefficient  (where  we  assume  damping  in  a 
span  is  linear,  viscous  damping),  and  m  is  the  beam  mass  per  unit  length.  yk  is 
displacement  of  the  k-th  beam  where  the  vehicle  is  traveling.  Fk(x,t)  is  the 
exciting  force  of  the  k-th  beam  due  to  the  multicar,  multiload  vehicle  acting  on  the 
beam, 


Fk(x,t)  -  £fki(t)5(xkj  -vt), 
k  !=  1 


(14) 


fk.(t)-  cp(ypjj  yki)  +  kp(ypjj  Ykj )] » 


where  yPij  is  the  displacement  of  primary  suspension  of  i-th  car  and  j-th  magnet 
on  the  k-th  beam,  ykj  is  the  displacement  of  k-th  beam  on  the  point  kj 
corresponding  to  the  displacement  ypjj,  and  kn  is  the  total  number  of  forces 
applied  to  the  k-th  beam  by  the  vehicle. 
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For  simply  supported  beams,  the  boundary  conditions  of  the  k-th  beam  are 


yk(t,o)=i^M»  =  o, 

OX 

r  \  d2y(t,L)  „ 

yk^D-  ^2 

ox 


(16) 


If  there  is  a  double-span  beam  (total  length  is  2L),  the  slope  and  bending  moment 
at  an  interior  simple  support  must  be  continuous  (Cai  et  al.  1992a);  thus 

yk<t,x)UL-  =  yk(t,x)|x^L+ =o, 


dyk(t,x) 


dx 


_  dyw(t,x) 


x>L- 


dx 


(17) 


x-»L+ 


d2yk(t,x) 


dx 


d  yk(t,x) 


x — — 


dx 


x-»L+ 


and  there  are 


yk(t,2L)=32^(t’22L)  =  0. 


dx^ 


(18) 


The  initial  conditions  are 


yk(x,0)  =  ^°U. 


(19) 


In  the  modal  analysis  method,  displacement  of  the  beam  is  expressed  as 


yk(x,t)  =  XQkn(t)<Pn(x), 

n=l 


(20) 


where  qkn(t)  are  time-varying  modal  amplitudes  and  <pn(x)  are  modal  shape 
functions  that  are  orthogonal  over  the  beam  length  0  <  x  <  L.  For  a  single-span 
beam, 


<pn(x)  =  V2sin|^^-j  =  V2sin^Xn-jr-j,  n  =  1,2,3,...; 


(21) 


2.4  Numerical  Simulations 

Numerical  simulations  of  dynamic  interactions  of  vehicle/guideway  systems, 
schematically  shown  in  Figs.  1  and  2,  were  carried  out  on  the  basis  of  the 
governing  equations  described  in  Sections  2.2  and  2.3  for  the  vehicle  and 
guideway.  Because  of  the  coupled  dynamic  interaction  between  the  vehicle  and 
guideway  (as  indicated  in  Eq.  10  where  guideway  deflections  are  input  to  the 
vehicle,  and  in  Eq.  15  where  vehicle  static  weight  and  acceleration  forces  are 
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excitations  to  the  guideway),  an  iterated  method  is  required  in  numerical 
simulations  to  calculate  dynamic  response  of  both  vehicle  and  guideway,  when 
the  fourth-order  Runge-Kutta  method  is  applied  in  the  simulations.  For  maglev 
vehicles  restricted  to  vertical  accelerations  of  less  than  0.05  g,  the  inertia  force  is 
much  lower  than  the  static  load,  and  dynamic  coupling  will  be  weak  (Richardson 
and  Wormley  1974).  In  this  case,  the  iteration  is  not  needed.  Because  the 
integrating  time-step  is  small  enough,  deflections  of  guideway  spans  in  the 
previous  time-step  can  be  used  as  input  to  the  vehicle,  and  dynamic  responses  of  a 
vehicle  can  then  be  calculated  and  the  results  used  to  calculate  guideway 
response  at  the  current  time-step.  This  calculating  sequence  proved  efficient 
when  coupling  between  the  vehicle  and  guideway  is  weak  or  when  vehicle  speed  is 
below  certain  values  (Cai  et  al.  1992a,  1992c). 

In  the  modal  analysis  method,  the  infinite  sum  in  Eq.  20  is  truncated  and  a 
finite  number  of  modes  is  used  to  represent  guideway  motion.  The  number  of 
modes,  n,  required  for  a  given  level  of  accuracy  depends  on  the  frequency  content 
of  the  guideway  forcing  function  (Eqs.  14  and  15),  the  traverse  speed  of  this  forcing 
function,  and  the  beam  properties.  For  maglev  systems,  the  number  of  modes  is 
related  to  vehicle  speed  or  crossing  frequency  (Cai  et  al.  1992a;  Chiu  et  al.  1971; 
Richardson  and  Wormley  1974).  Five  modes  are  sufficient  for  maglev  guideways 
with  vehicle  speeds  under  500  km/hr,  in  accordance  with  calculated  results  (Cai 
et  al.  1992a). 

The  focus  of  our  study  is  the  steady-state  or  repetitive  condition  of  guideway 
deflections  and  vehicle  heave  accelerations  for  the  vehicles  with  a  vertical  motion. 
The  steady  state  exists  after  a  vehicle  with  a  given  arbitrary  set  of  initial 
conditions  has  traversed  a  sufficient  number  of  spans  in  which  the  state  of  the 
vehicle  entering  a  span  is  identical  to  its  state  when  leaving  the  span  or,  in  fact, 
entering  the  next  span.  For  a  vehicle  starting  under  zero  initial  conditions,  the 
number  of  spans  a  vehicle  must  cross  to  reach  a  steady-state  condition  depends  on 
the  number  of  modes  and  traveling-speed  ratio  of  the  vehicle.  The  maximum 
number  of  spans  a  vehicle  must  cross  to  reach  a  steady  state  is  less  than  100,  in 
accordance  with  calculated  results  (Cai  et  al.  1992a). 

In  addition,  to  avoid  numerical  overflow  ir  simulation  from  the  effect  of 
equally  spaced  multiple  cars  that  may  excite  a  very  large  resonance  of  guideway 
deflections  (Richardson  and  Wormley  1974),  the  length  of  each  car  should  not  be 
the  same  as  the  length  of  the  beam  spar 
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2.5  Dynamic  Analyses  of  Vehicle/Guideway  interactions 

Table  1  shows  the  parameters  of  vehicle  and  guideway  we  used  in  our 
simulation  for  the  maglev  systems  shown  in  Figs  1  and  2. 

Two-Degree-of-Freedom  Vehicle.  If  we  set  M  =  1  and  N  =  1  in  Eqs  ]  1  and 
12,  the  vehicle  appears  to  be  a  two-degree-of-freedom  model  that  provides  a 
relatively  simple  explanation  of  the  dynamic  behavior  of  vehicle/guideway 
systems  (Cai  et  al.  1992a,  1992c). 

Figures  3  and  4  show  the  time  histories  of  steady-state  guideway 
displacement  ratio  Yg  (=  yg/y,n )  and  vehicle  acceleration  ratios  Yp  (=  yp/ym>  and 
Ys  (=  ys/ym)  for  both  primary  and  secondary  suspensions  on  the  single-span  beam 
vs.  the  location  of  the  vehicle  on  span  4v  1=  vt/L),  for  two  vehicle-traveling-speed 
ratios,  i.e.,  v/vc  =  0.25  in  Fig.  3  and  v/vc  =  0.5  in  Fig.  4,  where  vc  is  critical  speed 
and  is  equal  to  2fiL  (f[  -  coi/2tc).  ym  is  the  midspan  displacement  associated  with 
the  fundamental  mode  when  a  concentrated  static  load  (m.)  +  mu)g  is  placed  at  the 
midspan,  i.e.,  ym  =  2(mp  +  ms)g/mLu>J\  The  nondimensional  parameters  in  this 
model  are  chosen  as  {mp  +  ms)/(mL)  =  0.5;  mp/ms  =  0.1;  wi/ojp  =  0.3;  (os/o>p  =  0.25; 

=  2%;  £p  =  10%;  Cs  =  25%.  The  results  show  that  dynamic  interaction  between 

Table  1.  Parameters  for  dynamic  interaction  analysis  of  maglev  systems 
Vehicle 


Vehicle  length  l 

25.0  m 

Magnet  mass  mp 

1016  kg 

Car  body  mass  ms 

45700  kg 

Primary  damping  cp 

3.45  x  104  N-s/m 

Secondary  damping  cs 

2. 15  x  104  N-s/m 

Primary  stiffness  kp 

1.45  x  104  N/m 

Secondary  stiffness  ks 

2.26  x  104  N/m 

Intercar  vertical  stiffness  kv 

2.26  x  104  N/m 

Intercar  vertical  dampings  cv 

0.0 

Length  of  span  L 

25.0  m 

Bending  rigidity  El 

7.16  x  109  N-m2 

Mass  per  unit  length  m 

1.82  x  103  kg/m 

Damping  Cn 

3% 

Accat.  of  secondary  suspension,  Accal.  of  primary  suspension,  V  Olsplacamant  ratio  of  guldaway, 
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v/vc  .  0.25 
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Location  of  vehicle, 


Fig.  3.  Time  histories  of  steady-state  guideway 
displacement  ratio  and  vehicle 
acceleration  with  speed  ratio  v/vc  =  0.25 
for  a  two-degree  -  of -freedom  vehicle 


Accal.  ol  secondary  suspension,  ¥#  Accsl.  of  primary  suapansion,  V  Displacement  ratio  of  guideway. 


V/Ve  .  0.5 
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•a —  uncoupling 
-a—  coupling 


Location  of  vehicle, 


Location  of  vahlcla, 


Location  of  vahlcla,  4v 

Fig.  4.  Time  histories  of  steady-state  guideway 
displacement  ratio  and  vehicle 
acceleration  with  speed  ratio  v/vc  =  0.50 
for  a  two-degree-of- freedom  vehicle 
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vehicle  and  guideway  (comparing  uncoupled  and  coupled  curves)  has  relatively 
little  influence  on  the  secondary  suspension  at  the  given  parameter  values;  the 
effect  on  guideway  displacement  is  smaller  for  v/vc  =  0.25  than  for  v/vc  =  0.5,  but 
the  effect  on  the  acceleration  of  the  primary  suspension  is  greater  for  v/vc  =  0.25 
than  for  v/vc  =  0.5.  Also,  we  found  that  the  location  of  vehicle  at  maximum 
displacement  varies  with  the  vehicle-traveling-speed  ratio  v/vc.  In  Fig.  3,  v/vc  is 
0.25  and  is  about  0.4,  while  in  Fig.  4,  v/vc  is  0.5  and  is  about  0.65. 

Figure  5  shows  the  maximum  guideway  displacement  ratio  Yg  (=  yg/ym)  and 
maximum  vehicle  acceleration  ratios  Yp  (=yp/ym)  and  Ys(-ys/ym)  for  both 
primary  and  secondary  suspensions  as  a  function  of  vehicle-traveling-speed  ratio 
v/vc  on  both  single-  and  double-span  guideways.  For  a  single-span  beam,  the  peak 
of  maximum  displacement  ratio  is  about  1.7  when  v/vc  is  about  0.6  for  an 
uncoupling  model,  while  for  a  coupling  model  the  peak  of  maximum 
displacement  declines  to  1.5  when  v/vc  is  about  0.4.  When  v/vc  is  below  0.4,  the 
coupled  and  uncoupled  modes  remain  in  good  agreement.  Therefore,  for  small 
values  of  v/vc,  an  uncoupling  model  may  be  sufficient  to  simulate  dynamics  of 
vehicle/guideway  systems,  i.e.,  dynamic  motions  of  the  vehicle  and  guideway  can 
be  decoupled,  the  vehicle  may  be  considered  a  moving  force  on  the  guideway,  and 
guideway  deflection  is  then  used  as  a  known  displacement  input  into  the 
suspensions.  For  the  double-span  guideway,  when  the  vehicle-traveling-speed 
ratio  v/vc  <  0.5,  the  maximum  displacement  ratios  at  both  %  =  0.5  and  £,  =  1.5  (£  = 
x/L)  are  much  smaller  than  those  of  the  single-span  guideway.  The  differences 
between  uncoupling  and  coupling  models  for  the  double-span  guideway  are 
smaller  than  those  of  the  single-span  guideway.  From  comparisons  of  vehicle 
accelerations,  the  amplitudes  of  maximum  accelerations  of  both  primary  and 
secondary  suspensions  for  the  double-span  guideway  are  lower  than  those  of  the 
single-span  guideway.  Because  v/vc  in  maglev  systems  is  expected  to  be  no  larger 
than  0.5  (Sinha  1987),  and  without  considering  other  factors,  a  two-span  beam 
appears  to  be  more  efficient  in  achieving  better  ride  quality. 

More  detailed  parameter  analyses  for  the  two-degree-of-freedom  vehicle 
model  can  be  found  in  Cai  et  al.  1992a. 

Effects  of  Distributed  Loads.  In  a  dynamic  analysis  of  vehicle/guideway 
interactions,  an  understanding  of  the  effects  of  distributed  loads  is  essential.  In  a 
single-car  vehicle  (system  parameters  are  given  in  Table  1)  as  shown  in  Fig.  2,  for 
any  given  span  configuration,  span  deflections  decrease  as  the  number  of 
magnets  is  increased  and  total  force  is  held  constant.  These  effects  exist  when  the 
vehicle  travels  at  certain  speeds.  Figure  6  shows  the  midspan  deflections  of  a 
single-span  beam  when  a  single-car  vehicle,  which  has  one,  two,  four,  and  eight 
magnets  attached,  travels  at  100  m/s.  Figure  7  shows  the  maximum  midspan 
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v/vc 


v/vc 


Fig.  5.  Effects  of  single-span  and  double-span 
beams  on  maximum  guideway  displace¬ 
ment  and  vehicle  acceleration  ratios  for  a 
two-degree-of- freedom  vehicle 
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Fig.  6.  Midspan  displacement  of  guideway 
when  a  single  car  with  various 
magnets  travels  along  the  guideway  at 
100  m/s 


Traveling  Velocity,  m/s 

Fig.  7.  Maximum  midspan  displacements  of 
guideway  when  a  single  car  with 
various  magnets  travels  along  the 
guideway  at  various  speeds 

deflections  as  a  function  of  vehicle  traveling  velocity.  Apparently,  the  one-magnet 
case,  which  represents  a  two-degree-of-freedom  vehicle  with  a  concentrated  load, 
causes  the  largest  beam  deflection.  The  responses  of  four  magnets  and  eight 
magnets  have  almost  the  same  order  deflections  when  the  traveling  velocity  is 
greater  than  50  m/s. 
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Figure  8  shows  car-body  vertical  accelerations  with  one,  two,  four,  and  eight 
magnets  attached  to  the  vehicle.  Note  that  distributed  loads  affect  absolute 
acceleration  levels.  We  have  conducted  many  calculations  to  increase  the  number 
of  magnets  above  eight,  but  the  results  are  the  same  as  those  with  eight  magnets. 
Therefore,  we  believe  that  eight  magnets  are  sufficient  for  accurate  modeling  of 
distributed  loads  in  the  simulation.  Thus,  we  always  use  eight  magnets  for  each 
car  when  simulating  dynamics  of  a  multiple-car  vehicle  (described  later). 

Figures  9(a)  and  (b)  show  the  steady-state  accelerations  of  car  body  and 
magnets,  respectively,  when  a  single  car  with  eight  magnets  travels  along  the 
guideway  at  100  m/s. 

Dynamics  of  Multicar  Vehicle.  Simulations  on  dynamics  of  a  multicar 
vehicle  are  completed  by  using  the  model  given  in  Fig.  1. 

Figure  10  shows  midspan  beam  deflections  when  multicar  vehicles  (1,  2,  3, 
and  4  cars)  travel  at  100  m/s.  No  matter  how  many  cars  are  included  in  the 
vehicle,  the  maximum  beam  deflection  remains  the  same.  But  the  duration  of 
deflections  increases  as  car  number  increases.  Figure  11  shows  the  maximum 
displacements  of  the  guideway  midspan  when  the  multicar  vehicle  travels  at 
various  speeds.  Again,  results  for  1,  2,  3,  and  4  cars  are  the  same.  As  in  previous 
studies  on  the  concentrated-load  single-car  vehicle,  maximum  guideway 
displacements  tend  to  increase  as  vehicle  speed  increases. 


Fig.  8.  Accelerations  of  car  body  when  a  single 
car  with  various  magnets  travels  along 
the  guideway  at  100  mis 
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Fig.  9(a).  Acceleration  of  car  body  when  a 
single  car  with  eight  magnets 
travels  along  the  guideway  at 
100  mis 
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Fig.  9(b).  Accelerations  of  first  and  eighth 
magnets  when  a  single  car  with 
eight  magnets  travels  along  the 
guideway  at  100  m/s 

Figure  12  shows  car  body  accelerations  for  vehicles  with  various  cars  when 
traveling  speed  is  100  m/s.  We  note  that  the  single-car  vehicle  has  the  largest 
peak-to-peak  acceleration  and  that  the  multicar-vehicle  peak-to-peak  acceleration 
decreases,  which  indicates  that  intercar  restraints  affect  vehicle  motions  and  that 
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Fig.  10.  Midspan  displacements  of  guideway 
for  multicar  vehicles  with  eight 
magnets  on  each  car  traveling  along 
the  guideway  at  100  mis 


Traveling  Velocity,  m/s 

Fig.  11.  Maximum  midspan  displacements  of 
guideway  when  multicar  vehicles 
with  eight  magnets  on  each  car  travel 
along  the  guideway  at  various  speeds 

the  multicar  vehicles  may  have  better  ride  comfort.  Figure  13  gives  the  peak-to- 
peak  car  body  accelerations  for  vehicles  at  various  speeds;  the  results  are  the 
same  as  in  Fig.  12.  The  peak-to-peak  accelerations  for  a  single  car  are  much 
larger  than  those  of  a  multicar  vehicle. 
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Fig.  12.  Accelerations  of  car  body  for  multicar 
vehicles  with  eight  magnets  on  each 
car  traveling  along  the  guideway  at 
100  m/s 


Traveling  Velocity,  m/s 


Fig.  13.  Peak-to-peak  accelerations  of  car  body 
when  multicar  vehicles  with  eight 
magnets  on  each  car  travel  along  the 
guideway  at  various  speeds 

To  determine  effects  of  intercar  constraints  on  vehicle  dynamics,  further 
calculations  are  completed  for  different  intercar  stiffness,  and  the  results  for  a 
two-car  vehicle  are  shown  in  Fig.  14.  As  stiffness  decreases,  peak-to-peak 
acceleration  increases. 
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Fig.  14.  Accelerations  of  car  body  when  a  two- 
car  vehicle  with  eight  magnets  on 
each  car  travels  along  the  guideway  at 
100  m/s  and  when  intercar  stiffness 
varies 

Ride  Comfort  of  Muiticar  Vehicle.  Figure  15  shows  a  comparison  of  the 
Urban  Tracked  Air  Cushion  Vehicle  (UTACV)  ride  comfort  specification 
(ranging  from  0-10  Hz)  for  multicar  vehicles  traveling  at  100  m/s.  With  the 
parameters  from  Table  1,  power  spectral  densities  (PSDs)  of  multicar  vehicle 
accelerations  satisfy  the  ride  comfort  criterion.  It  appears  that  the  vehicle  with 
those  parameters  can  provide  an  acceptable  ride.  From  Fig.  15,  we  also  note  that 
at  the  fundamental  frequency,  the  PSDs  of  acceleration  decrease  as  car  number 
increases;  however,  at  higher  harmonic  frequencies,  this  tendency  is  not  so  clear. 

Figures  16  and  17  show  PSDs  of  acceleration  of  a  single-car  and  a  two-car 
vehicle,  respectively,  traveling  at  various  speeds;  the  harmonic  frequencies  vary 
with  traveling  speed. 


3  Control  Designs  for  Maglev  Systems _ 

3.1  Background  of  Control  Law  Designs 

Suspension  systems  are  dominant  in  determining  the  basic  dynamic  and 
vibrational  behavior  of  the  maglev  vehicle.  Primary  suspensions,  which  consist  of 
magnets,  provide  basic  support  and  guidance  along  the  guideway  and  should 
respond  quickly  to  guideway  disturbances  and  perturbations.  Secondary 
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Fig.  15.  PSD  of  car  body  accelerations  when 
multicar  vehicles  with  eight  magnets 
on  each  car  travel  along  the  guideway 
at  100  m/s 


Fig.  16.  PSD  of  car  body  accelerations  when  a 
single-car  vehicle  with  eight  magnets 
travels  along  the  guideway  at  various 
speeds 

suspensions  support  the  vehicle  bodies  and  provide  acceptable  ride  comfort  to  the 
passengers. 

To  achieve  a  quick  response  and  a  high-quality  ride  over  a  less-expensive 
guideway,  control  designs  must  be  exploited  in  suspension  systems.  Moreover, 
with  the  assistance  of  suspension  controls,  a  rougher  guideway  surface  could  be 
used  and  overall  investment  cost  of  the  guideway  could  be  reduced. 
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Fig.  1 7.  PSD  of  car  body  accelerations  when  a 
two-car  vehicle  with  eight  magnets  on 
each  car  travels  along  the  guideway  at 
various  speeds 

Primary  suspension  control  strategies  include  two  basic  principles:  position 
control  and  air-gap  control.  The  position  control  maintains  the  vehicle  in  a  stable 
levitation  and  guidance  attitude  against  the  various  guideway  irregularities  that 
are  the  dominant  excitation  force  to  the  suspension  systems  and  may  be 
represented  by  a  stationary  stochastic  signal.  This  control  guides  the  vehicle  on 
the  guideway  curve  purely  on  the  basis  of  alignment  and  compensates  for  all 
deviations  of  the  guideway  from  this  ideal  line.  Primary  suspension  controls 
must  also  respond  quickly  to  small  air-gap  changes  and  overcome  the  effects  of 
high-frequency  perturbations.  The  air-gap  control  aims  at  following  the  actual 
path  of  the  track,  including  its  deviations  from  the  line  of  alignment  caused  by 
constructional  features,  to  maintain  a  constant  air-gap. 

Secondary  suspension  control  systems  should  provide  ride  comfort  to  the 
passengers  over  guideway  sections  that  are  irregular  or  when  the  vehicle  is 
operating  in  gusty  winds. 

Even  though  suspension  control  designs  are  important  to  maglev  systems, 
only  a  few  studies  on  the  control  design  of  maglev  suspensions  can  be  found  in  the 
literature.  Katz  et  al.  (1974)  studied  a  linearized  control  system  for  magnets  with 
a  proportional  feedback  in  primary  suspension.  Kortum  and  Utzt  (1984)  used 
quadratic  synthesis  in  the  control  design  of  a  vertical  maglev  vehicle  model. 
Kortum  et  al.  (1988)  and  Faye  et  al.  (1989)  studied  dynamic  modeling  of  maglev 
vehicles  for  control  design  and  performance  evaluations.  They  applied  two 
control  strategies,  i.e.,  feedback  controls  with  pole  assignment  and  with  the 
Riccati  (quadratic)  design. 
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Existing  maglev  suspension  systems  contain  passive  control  elements  such 
as  dampers  and  springs  (Bohn  and  Alscher  1986;  Bohn  and  Steinmetz  1985), 
which  are  passive  in  the  sense  that  no  power  source  is  required,  i.e.,  the 
vibration-control  elements  only  store  or  dissipate  the  energy  associated  with 
vibration.  Such  passive  suspension  systems  may  be  insufficient  to  improve 
suspension  performance  unless  radical  changes  are  made  to  the  basic  principle 
of  suspension  operation  (Faye  et  al.  1989;  Kortum  et  al.  1988). 

The  alternative  proposed  here  is  that  of  an  active  suspension  system,  which 
has  been  investigated  by  many  researchers  over  the  last  30  or  more  years  on 
conventional  ground  vehicles  due  to  its  potential  to  improve  vehicle  performance 
(Bernard  et  al.  1987;  Chalasani  1987;  Elmadany  1990;  Dukkipati  et  al.  1992;  Goodal 
and  Kortum  1983;  Karnopp  et  al.  1974;  Karnopp  and  Margolis  1984;  McCormac  et 
al.  1992). 

An  active  suspension  system  offers  a  solution  to  the  conflicting  requirements 
of  a  constant-parameter  suspension  design  and  allows  continuous  or  discrete 
variation  in  effective  spring  constants  and  damping  coefficients.  An  active 
suspension  may  be  adapted  to  specific  vehicle/guideway  conditions  by  control 
devices  that  can  adjust  the  dynamics  of  hydraulic  or  pneumatic  actuators.  In 
such  a  system,  the  passive  spring  and  damper  of  a  conventional  suspension  are 
replaced  by  force  actuators  that  continuously  supply  and  modulate  the  flow  of 
energy  by  generating  forces  on  the  basis  of  some  control  law  to  achieve  the 
required  performance  (Bernard  et  al.  1987).  It  is  apparent  that  an  active  system 
exhibits  significantly  low  PSD  values.  In  other  words,  the  active  system  does  not 
amplify  accelerations  as  one  would  expect  with  a  passive  system.  Also,  in  the 
active-system  concept,  algorithms  can  be  designed  (in  software  rather  than 
hardware)  to  provide  the  same  function  to  achieve  suspension  performance 
(McCormac  et  al.  1992).  In  general,  however,  active  systems  are  more  costly, 
more  complex,  and  often  less  reliable  than  passive  systems  (Bernard  et  al.  1987; 
Karnopp  et  al.  1974). 

Semiactive  suspensions,  which  have  been  under  development  in  recent  years 
(Alonoly  and  Sankar  1987,  1988;  Hrovat  et  al.  1988;  Karnopp  et  al.  1974),  present  a 
compromise  between  active  and  passive  suspensions.  Semiactive  suspensions 
can  achieve  performance  close  to  that  of  active  suspensions  with  much  lower  cost 
and  complexity.  These  systems,  which  require  relatively  little  external  power, 
use  an  active  damper  in  parallel  with  a  passive  spring.  Desired  forces  in  the 
damper  are  generated  by  modulating  fluid-flow  orifices  via  a  control  scheme 
involving  feedback  variables. 

Although  semiactive  dampers  based  upon  the  electro-hydraulic  servovalve  do 
not  consume  as  much  power  as  fully  active  systems,  they  do  require  a  number  of 


mechanical  components,  each  manufactured  to  fine  tolerances.  One  way  to 
simplify  the  construction  of  a  semiactive  device  is  to  obtain  continuous  control  of 
damping  by  exploiting  "smart"  materials.  Among  these  different  materials,  the 
electro-rheological  fluids  (ER  fluids),  which  have  variable-shear  characteristics, 
can  undergo  significant  instantaneous  reversible  changes  in  damping 
characteristics  when  subjected  to  electrostatic  potentials.  This  makes  their 
application  to  real-time  semiactive  vibration  control  very  attractive.  ER  fluids 
were  discovered  by  Willis  M.  Winslow  in  1939  and  have  attracted  intense  interest 
in  recent  years  (Stanway  et  al.  1989;  Wong  et  al.  1992;  Wu  et  al  1990).  Because  the 
apparent  viscosity  of  ER  fluids  can  be  conveniently  controlled  by  an  applied 
electric  field,  permitting  the  fluids  to  change  from  normal  to  viscous  and  back 
again  in  less  than  a  millisecond,  the  damping  effect  can  be  altered  even  at  high 
speeds. 

It  is  thought  that  ER  fluids  are  particularly  suited  for  continuously 
adjustable  dampers  in  semiactive  suspension  systems  for  ground  vehicles.  In 
essence,  this  type  of  damper  behaves  like  a  conventional  viscous  damper  with  an 
adjustable  orifice,  but  is  controlled  by  activating  the  electrodes  and  thus  needs  no 
mechanical  moving  parts.  In  addition,  such  devices  could  be  made  without  the 
need  for  precision  machining,  special  fits,  or  exotic  materials.  It  should  be 
possible  to  reduce  component  costs,  in  relation  to  active  system,  by  employing 
continuously  controllable  damping  available  with  ER  fluids.  The  use  of  ER  fluids 
as  a  direct  interface  between  mechanical  suspension  components  and  control 
electronics  offers  an  elegant  solution  to  the  problem  of  implementing  semiactive 
control.  A  semiactive  suspension  system  with  ER  dampers  employing  the 
continuous  control  strategy  has  the  potential  to  provide  ground  vehicles  with 
significantly  improved  ride  comfort  over  that  of  conventional  passive  suspensions 
(Wong  et  al.  1992). 


3.2  Modeling  of  Maglev  Vehicle  Suspensions  for  Control  Designs 

To  investigate  the  improvement  of  the  dynamic  response  and  ride  comfort  of 
maglev  systems,  different  control  designs  (active  and  semiactive)  are  examined  in 
this  study. 

For  most  control  law  synthesis,  it  is  desirable  to  work  with  linear  dynamic 
models  of  low  order  A  low-order  maglev  vehicle  model,  which  may  be  selected  as 
a  two-degree-of-freedom  quarter-vehicle  model  representing  primary  and 
secondary  suspensions,  is  necessary  in  control  design  to  formulate  a  low'-order 
controller.  Again,  a  vertical  model  is  sufficient  to  perform  control-law  design  for 
maglev  systems  when  considering  vehicle  dynamics  and  ride  comfort  to 
passengers. 
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For  these  reasons,  a  one-dimensional  vehicle  model  with  two  degrees  of 
freedom  (Fig.  18)  and  consisting  of  two  lumped  masses  mp  and  ms,  two  linear 
springs  kp  and  ks,  and  two  viscous  dampings  cp  and  cs,  representing  primary  and 
secondary  suspensions,  respectively,  is  used  in  the  control  synthesis  of  maglev 
systems.  This  model  is  simple  to  analyze  and  is  used  mainly  to  gain  a  basic 
understanding  of  the  suspension  concept.  Such  a  simplified  model  can  be 
sufficient  to  predict  the  behavior  of  the  system.  In  addition,  with  this  model,  the 
coupling  effects  of  primary  and  secondary  suspensions  will  not  be  ignored;  this  is 
very  important  for  suspension  control  design  because  coupling  of  both 
suspensions  may  result  in  interaction  forces  or  kinematic  constraints  and  can 
thus  affect  ride  comfort, 

A  mechanical  control  method  is  favored  in  our  study  because  it  is  able  to 
change  conjugate  eigenvalues  of  the  system,  while  controlled  coils  or  controlled 
magnet  current  can  only  change  real  poles  of  the  system  (Faye  et  al.  1989;  Kortum 
etal.  1988). 

The  passive  parameters  of  the  German  Transrapid  Maglev  System  TR06, 
(summarized  in  Table  2)  are  utilized  for  analysis  in  this  study  because  no  other 
Transrapid  data  are  available  in  the  literature  (Bohn  and  Steinmetz  1985;  Bohn 
and  Alscher  1986). 


Fig.  18.  One-dimensional  two-degree-of- 
freedom  vehicle  model  with 
primary  and  secondary  suspen¬ 
sions  for  maglev  systems 
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Table  2.  Vehicle  vertical  dynamics  model  parameters  of  German 
Transrapid  Maglev  System  TR06 


Primary  suspension  mass  mp 

3.20  x  104  kg 

Secondary  suspension  mass  ms 

2.92  x  104  kg 

Primary  suspension  damping  cp 

1.13  x  106  N-s/m 

Secondary  suspension  damping  cs 

8.80  x  104  N-s/m 

Primary  suspension  stiffness  kp 

6.18  x  107  N/m 

Secondary  suspension  stiffness  ks 

7.37  x  105  N/m 

For  this  model,  the  equations  of  motion  when  the  vehicle  is  at  an  equilibrium 
position  are  then 

mPyP + cP(yP  -  yi) + kP(yP  -y0~  cs(ys  -  yP)  -  Mys  -  yP>  =  °>  <27> 

msys+cs(ys-yP)  +  ks(ys-yp)  =  0,  (28) 

from  Eqs.  27  and  28,  there  is 

mpyp  +  msys  +  cp  (yp  -  yi )  +  kp  (yp  -  y* )  =  0,  (29) 

where  yp  and  ys  are  positions  of  primary  and  secondary  suspensions,  and  yj  is 
disturbance  from  guideway  irregularity.  System  parameters,  i.e.,  masses  mp 
and  m,,  stiffnesses  kp  and  ks,  and  dampings  cp  and  cs  for  primary  and  secondary 
suspensions,  are  listed  in  Table  2. 

Table  3  shows  eigenvalues,  frequencies,  and  damping  ratios  for  both  primary 
and  secondary  suspensions  in  this  model. 


3.3  Dynamic  Response  of  Maglev  Suspensions 

Taking  the  Laplace  transform  of  the  differential  equations  of  Eqs.  28  and  29, 
and  assuming  all  initial  conditions  to  be  zero,  yields 

|mps2  +cps  +  kpjYp  +mss2Ys  -(cps  +  kp)Yj  =  0  (30) 


and 
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Table  3.  Eigenvalues ,  frequencies ,  and  damping  ratios  in  1-D  vehicle 
model  with  two  degrees  of  freedom 


Eigenvalue 
o±  jco 

Frequency 
f  (Hz) 

Damping 

z 

TR06 

-1.479  ±4.779j 

0.80 

0.30 

-19.06  ±  39.8  lj 

7.03 

0.43 

Passive  Adjustments 

-3.024  ±  4.027j 

0.80 

0.60 

-22.02  ±  37.93j 

6.98 

0.50 

Primary  Feedback 

-1.479  ±  4.779j 

0.80 

0.30 

-23.89 

— 

— 

-65.71  ±84.55j 

17.05 

0.61 

Primary  and 

-3.014  ±  4.020j 

0.80 

0.60 

Secondary 

-23.89 

— 

— 

Feedback 

-65.71  ±  84.55j 

17.05 

0.61 

(mss2  +  css  +  ks)Ys  -(css  +  ks)Yp  =0,  (31) 

where  s  is  the  Laplace  transform  operator,  and  Yp,  Ys,  and  Yj  are  the  Laplace 
transforms  of  yp,  ys,  and  yj. 

From  Eqs.  30  and  31,  transfer  functions  for  primary  and  secondary 
suspensions  when  considering  guideway  disturbance  Yj  as  an  input  variable  can 
be  as  described  below: 


Gp(s)  = 


Yp(s) 

Yj(s) 


(cps  +  kp)(mss2  +  css  +  ks) 

(mps2  +  cps  +  kpj|mss2  +  css  +  ks )  +  mss2(css  +  ks ) 

(^pkg  cskp  js  +  kpkg 

mpmss4  +  [mpcs  +  ms(cp  +  cs)]s3  +  [mpks  +  ms(kp  +  ks)  +  cpcs]s2  +  (cpks  +  cskp)s  +  kpks 


(32) 
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Gs(s) 


Ys(s) 

Yj(s) 


= _ (cps  +  kp)(c8s  +  kg) _ 

(mps2  +  cps  +  kp)(mss2  +  c8s  +  k8)  +  mas2(css  +  k8) 

_  CpCss  +(cpks  +cskp)s+kpkg 

mpm6s4  +[mpcs  +  ms(cp  +c8)]s3  +(mpk8  +  m8(kp  +  k8)  +  cpc8Js2  +{cpk8  +c8kp)s+  kpkB 

(33) 

From  Eqs.  30  to  33,  by  connecting  signals  properly,  we  can  construct  a  block 
diagram  for  the  system,  as  shown  in  Fig.  19. 

In  many  practical  cases,  the  desired  performance  characteristics  of  control 
systems  are  specified  in  terms  of  time-domain  quantities  or  of  transient  response 
to  a  unit-step  input  because  such  input  is  easy  to  generate  and  is  sufficiently 
drastic.  The  transient  response  of  a  system  to  a  unit-step  input  depends  on  the 
initial  conditions.  For  convenience  in  comparing  transient  responses,  it  is 
common  practice  to  use  the  standard  initial  condition  that  the  system  is  at  rest 
initially,  with  output  and  all  time  derivatives  thereof  equal  to  zero. 

From  the  given  transfer  function  in  Eqs.  32  and  33,  the  corresponding 
differential  equation  can  be  obtained  as 

yp+aiyp  +  *2yp  +  a3yp  +  a4yp  =  bplyj  +  bp2yj  +  bp3yj  +  bp4yj  (34) 

ys+aiys  +  a2ys  +  a3ys  +  a4ys  =  bsiyi  +  tw*  +bs3yt  +  bs4yi,  (35) 


Fig.  19.  Block  diagram  for  two-degree-of -freedom  vehicle  model 
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where 


al 


mpcs  +  ms(cp+cs) 
mpms 


mpks  +  ms(kp  +  ks)  +  cpcs 

a2  = - 

mpms 


cpks  +  cskp 

a3  _ - 

mpms 


kpkg 
a4  =  ~ 


mpms 


(36) 


mscp  cp 
mpms  mp 


mskp  +  cpcs 
mpms 


cpks  +  csk 
mpms 


p 


kpks 

mpms 


(37) 


and 


We  define  state  variables  xi,  X2,  x3>  and  X4  as  follows: 


xi  =  yP 

x2=xi-piyi 

x3=x2-p2yi 

X4=x3-P3yi 

where 

Pi  =  bpl 
P2  =  bp2  -  aipi 
P3  =  bp3  -  aiP2  -  a2pi 
p4  =  bp4  -  aip3  -  a2p2  -  a3pi 
Note  that 


x4  =  -a4x !  -  a3x2  -  a2x3  -  a  xx4  +  p4  y j . 


(39) 


(40) 


(41) 


The  state  equation  and  output  can  be  obtained  as 
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(42) 


(43) 


Then  the  unit-step  response  of  primary  suspension  yp(t)  versus  time  t,  when  input 
variable  y\  reached  1,  can  be  calculated  from  Eqs.  42  and  43  by  the  fourth-order 
Runge-Kutta  method.  Similarly,  the  unit-step  response  of  the  secondary 
suspension  ys(t)  can  be  obtained  by  rewriting  state  variables  in  Eqs.  39  and  40. 

The  frequency  responses  of  the  vehicle  model,  by  which  we  mean  the  steady- 
state  responses  of  suspension  systems  to  a  sinusoidal  input,  can  be  obtained 
directly  from  the  transfer  function;  i.e.,  s  in  the  transfer  function  is  replaced  by 
jco,  where  co  is  angular  frequency.  Hence,  the  frequency  responses  of  Gp(jo))  and 
Gs(jco)  for  primary  and  secondary  suspension  with  respect  to  the  input  sinusoid 
can  be  obtained  from  Eqs.  32  and  33  by  substituting  jo)  for  s.  The  magnitude  ratios 
and  phase  shifts  of  Gp(jco)  and  Gs(jo))  can  be  presented  by  separate  plots  in  Bode 
diagrams,  that  is,  a  plot  of  the  logarithm  of  the  magnitude  of  a  sinusoidal  transfer 
function,  and  a  plot  of  the  phase  angle;  both  are  plotted  against  the  frequency  in 
logarithm  scale. 

With  the  frequency  transfer  function,  the  acceleration  output  PSD  of  the 
secondary  suspension,  which  is  a  measurement  of  vehicle  ride  comfort,  can  be 
derived.  The  surface  roughness  of  guideways  can  be  described  approximately  by 
the  PSD  of  surface  profile, 

S(co)  =  Av/o)2,  (44) 

where  v  is  vehicle  velocity,  and  A  is  roughness  amplitude  that  ranges  from  0.6  x 
10*6  to  20  x  10'6  m  (Gottzein  et  al.  1974;  Brock  1973).  Acceleration  output  PSD  of 
secondary  suspension  d>(o)  can  then  be  obtained  from  Katz  et  al.  (1974)  as 

O(o))  =  Av|j(oGs(jco)|2-  — , 

s 


(45) 
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where  g  is  acceleration  of  gravity,  and  the  factor  2k  is  needed  to  convert  from  a 
per-rad/s  to  a  per-Hertz  basis. 

With  the  data  given  in  Table  2,  dynamic  responses  of  the  system  (see  Figs.  18 
and  19)  are  then  calculated. 

Figures  20  and  21  show  transient  responses  of  primary  and  secondary 
suspensions  with  unit-step  input  of  guideway  perturbation.  Maximum 


Fig.  20.  Transient  response  yp  of  primary 
suspension  of  TR06  with  unit-step 
input  of  guideway  perturbation 


Fig.  21.  Transient  response  ys  of  secondary 
suspension  of  TR06  with  unit-step 
input  of  guide  way  perturbation 
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overshoots  of  Mp  for  primary  and  secondary  suspensions  are  31.6%  and  46.3%, 
respectively;  and  setting  times  of  ts  are  0.157  s  and  2.000  s  (5%  criterion  is  chosen); 
see  Tables  4  and  5.  Obviously,  these  dynamic  responses  are  not  desirable  for 
maglev  systems,  for  which  transient  responses  must  be  sufficiently  damped.  For 
a  desirable  transient  response  of  a  second-order  system,  the  damping  ratio  must 
be  between  0.4  and  0.8.  Based  on  data  in  Table  2,  the  damping  ratios  for  primary 
and  secondary  suspensions  are  0.43  and  0.30  (see  Table  3),  which  are  below  the 
desirable  value  and  yield  excessive  overshoots  in  transient  responses  (see  Figs.  20 
and  21,  and  Tables  4  and  5). 

Figures  22  and  23  show  frequency  responses  (Bode  diagrams)  for  primary 
and  secondary  suspensions;  these  indicate  particular  features  of  second-order 
systems  for  both  suspensions.  Responses  of  both  suspensions  have  excessive 
overshoots,  which  agrees  well  with  Figs.  20  and  21. 

Figure  24  shows  the  normalized  output  acceleration  PSD  of  vehicle  body  (ms), 
where  A  =  1  x  10*6  m  and  v  =  100  m/s.  From  Fig.  24,  the  first  peak  of  the  PSD  of 
vehicle  acceleration  satisfies  the  UTACV  ride  comfort  criterion;  however,  the 
second  peak  (which  corresponds  to  coupling  effects  of  the  primary  suspension) 
exceeds  this  criterion.  No  matter  what  other  features  are  used  in  the  secondary 
suspension,  oscillation  of  the  primary  suspension  still  affects  ride  comfort. 


Table  4.  Transient  measurement  of  unit-step  response  of  primary 
suspension  yp 


Mp(%) 

tp(s) 

ts(s) 

TR06 

31.6 

0.058 

0.157 

Passive  Adjustment 

26.8 

0.058 

0.104 

Primary  Feedback 

0.0 

— 

0.110 

Primary  and  Secondary  Feedback 

0.0 

— 

0.108 

Table  5.  Transient  measurement  of 

unit-step 

response  of  secondary 

suspension  ys 

Mp(%) 

tp(s) 

ts(s) 

TR06 

46.3 

0.520 

2.000 

Passive  Adjustment 

25.7 

0.456 

0.868 

Primary  Feedback 

44.8 

0.554 

2.020 

Primary  and  Secondary  Feedback 

12.3 

0.630 

0.964 
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Fig.  22.  Frequency  response  Gp(ja>)  of  primary 
suspension  of  TR06 


0.1  1.0  10.0 
Frequency,  Hz 

Fig.  23.  Frequency  response  Gs(jco)  of  secondary 
suspension  of  TR06 

3.4  Passive  Control  Design 

To  improve  dynamic  response  of  the  system,  a  passive  control  design  was 
examined.  The  passive  parameters  of  stiffness  and  dampings  were  utilized  in 
modulating  eigenvalues  of  the  system  to  provide  desirable  transient  response. 

As  shown  in  Table  3,  placing  two  pairs  of  conjugate  eigenvalues  of  the  system 
at  -3.024  ±4.027j  and  -22.02  ±37.93j  changes  the  damping  ratios  of  two 
suspensions  from  £p  =  0.43  and  £s  =  0.3  to  Cp  =  0.50  and  £s  =  0.60  while  system 
frequencies  remain  the  same. 
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Fig.  24.  PSD  of  vehicle  acceleration  of  TR06 
with  vehicle  speed  v  -  100  ml  sand 
guideway  roughness  amplitude 
A  -  106m 


Figures  25  and  26  show  the  transient  response  of  both  primary  and  secondary 
suspensions  with  unit-step  input  of  guideway  disturbance.  Maximum  overshoots 
of  Mp  for  primary  and  secondary  suspension  are  26.8%  and  25.7%,  respectively; 
and  setting  times  of  ts  are  0.104  s  and  0.868  s.  These  results  are  much  better  than 
that  of  TR06.  System  response  appears  to  be  faster  and  more  damped  than  before. 

Figures  27  and  28  are  frequency  responses  for  the  primary  and  secondary 
suspensions,  respectively,  and  show  the  same  good  results  as  in  Figs.  25  and  26. 

When  we  look  at  the  PSD  of  vehicle  acceleration  in  Fig.  29,  the  result  is  not 
satisfied  for  coupling  effects  of  the  primary  suspension  in  higher-frequency-range 
increase  and  ride  quality  becomes  worse.  This  is  because  when  we  change  poles 
of  transfer  function  of  the  system,  we  also  change  zeros  of  transfer  function  of  the 
system,  which  is  not  desirable  in  this  design.  It  is  verified  that  passive 
suspension  control  will  not  eliminate  acceleration  amplitude. 


3.5  Active  Control  of  Primary  Suspension 

As  mentioned  before,  primary  suspension  control  is  needed  to  maintain  the 
vehicle  in  a  stable  levitation  and  guidance  attitude  against  various  guideway 
irregularities  and  to  hold  a  constant  air  gap  between  magnets  and  guideway.  In 
addition,  the  primary  suspension  control  should  attenuate  coupling  effects  of 
primary  suspension  oscillation  during  acceleration  of  the  secondary  suspension, 
in  order  to  achieve  acceptable  ride  comfort.  With  a  passive  control  method, 
however,  these  goals  cannot  be  achieved  satisfactorily. 
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Fig.  25.  Transient  response  yp  of  primary 
suspension  with  unit-step  input  of 
guideway  perturbation  using 
passive  adjustments 


Fig.  26.  Transient  response  ys  of  secondary 
suspension  with  unit-step  input  of 
guideway  perturbation  using 
passive  adjustments 

An  active  primary  suspension  system  is  suitable  for  maglev  vehicles.  Such  a 
system  provides  continuous  or  discrete  variation  in  effective  spring  constants  and 
damping  coefficients,  according  to  some  control  law  that  may  be  designed  in 
software  rather  than  hardware.  Figure  30  shows  this  control  configuration  as  a 
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Fig.  27.  Frequency  response  Gp(jco)  of  primary 
suspension  using  passive  adjustments 


0.1  1.0  10.0  100.0 
Frequency,  Hz 


Fig.  28.  Frequency  response  Gs(jco)  of 
secondary  suspension  using  passive 
adjustments 

single- degree-of- freedom  model.  An  active  force  element  is  added  to  the  passive 
spring  and  damper  assembly  and  will  continuously  vary  the  force  acting  on  the 
primary  mass  to  change  the  response  characteristics  and  eliminate  acceleration 
amplitude.  In  this  approach,  the  force  element  can  be  realized  with  a  linear 
electrohydraulic  actuator  that  connects  magnet  and  bogie  in  the  primary 
suspension  (see  Fig.  31).  A  position  sensor  detects  the  air  gap  between  magnet 
and  guideway,  and  an  accelerometer,  mounted  on  the  bogie,  detects  bogie  motion. 
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Frequency,  Hz 


Fig.  29.  PSD  of  vehicle  acceleration  using 
passive  adjustments  with  vehicle 
speed  v  =  100  m/s  and  guideway 
roughness  amplitude  A  =  10  6  m 


Fig.  30.  Active  control  configuration  for  primary  suspension 

The  resulting  signals  are  processed  by  the  controller  according  to  designed 
control  law  in  software,  in  turn  causing  the  actuator  to  ensure  that  the  air  gap 
does  not  exceed  specific  tolerances  within  the  safety  margin  and  that  the 
acceleration  PSD  of  the  suspension  remains  as  low  as  possible  in  the  specified 
frequency  ranges  in  order  to  guarantee  good  ride  comfort.  In  practice,  when  the 
magnet  quickly  follows  guideway  deviations  caused  by  dynamic  and  static 
deflections  or  constructional  features,  the  magnet  oscillations  may  transfer  to  the 
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Accelerometer 
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Fig.  31.  Active  electro-hydraulic  system 

bogie.  Therefore,  an  optimal/adaptive  control  law  should  be  exploited  in  the 
controller,  and  some  trade-offs  may  be  needed  to  achieve  a  constant  air  gap  and  a 
low  acceleration  level  in  the  suspension  system. 

In  this  study,  an  active  feedback  control  path  is  applied  to  the  simplified 
vehicle  model  with  the  pole-assignment  approach,  which  provides  a  less 
complicated  control  model.  According  to  the  analysis  of  acceleration  PSD  in 
Fig.  24,  to  attenuate  the  coupling  effects  of  primary  suspension  oscillation  in  the 
high-frequency  range,  a  lag-lead  regulator  is  designed  in  the  inner  feedback  path 
for  the  primary  suspension  (see  Fig.  32). 

The  pole-assignment  method  involves  selecting  a  set  of  desired  pole  locations 
and  calculating  the  gains  required  to  achieve  such  a  pole  configuration  and  hence 
the  desired  system  response.  By  selecting  new  open-loop  zero  a  =  -5  (1/s),  and  pole 
b  =  -120  (1/s)  in  the  lag-lead  feedback  network,  the  root  locus  plot  of  the  primary 
suspension  is  given  in  Fig.  33.  For  feedback  gain  kp  =  2.7  x  108N,  the  closed-loop 
eigenvalues  are  given  in  Fig.  35  and  Table  3,  e.g.,  -23.89  and  -65.71  ±  84.55j;  this 
changes  the  damping  ratio  of  the  primary  suspension  to  0.61,  which  is  sufficient 
to  attenuate  the  oscillation,  and  the  frequency  to  17.05  Hz,  which  is  desirable 
because  the  higher  the  natural  frequency  of  the  primary  suspension  the  tighter 
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Fig.  32.  Block  diagram  for  two-degree-of -freedom  vehicle  model  of  maglev  system 
with  primary  and  secondary  suspension  feedback  controls 


a 

Fig.  33.  Root-locus  diagram  for  primary 
suspension  feedback  control  design 

the  tracking  will  be.  It  should  be  noted  that  with  mechanical  control  devices, 
such  as  electrohydraulic  servovalves,  the  complex  pair  of  suspension  eigenvalues 
can  be  easily  changed.  However,  the  controlled  coils  or  controlled  magnet  current 
can  change  only  the  real  poles  of  the  suspension  system. 

Figure  34  is  a  comparison  of  vehicle  acceleration  PSDs.  With  primary 
suspension  feedback  control,  the  second  peak  of  acceleration,  which  reflects  the 
coupling  effects  of  the  primary  suspension,  is  much  lower  and  meets  the  criterion 
with  parameter  selections  of  \p,  a,  and  b,  as  described  in  Fig.  32. 
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Fig.  34.  PSD  of  vehicle  acceleration  using 
active  feedback  controls  in  primary 
suspension  with  vehicle  speed  v  = 
100  m  I s  and  guideway  roughness 
amplitude  A  =  10  6  m 
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Fig.  35.  Transient  response  yp  of  primary  sus¬ 
pension  with  unit-step  input  of  guide¬ 
way  perturbation  using  active  feedback 
controls  in  primary  suspension 


Figures  35  and  36  compare  transient  responses  of  the  primary  and  secondary 
suspensions  with  unit-step  input  of  guideway  disturbance.  With  primary 
suspension  feedback  control,  dynamic  response  of  the  primary  suspension  is 
significantly  improved  (see  Fig.  35  and  Table  4)  and  the  setting  time  is  reduced  to 
ts  =  0.115  s,  but  it  does  not  affect  the  response  of  the  secondary  suspension  and  the 
setting  time  of  ts  remains  the  same  as  before  (see  Fig.  36  and  Table  5). 
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Fig.  36.  Transient  response  ys  of  secondary 
suspension  with  unit-step  input  of 
guideway  perturbation  using  active 
feedback  controls  in  primary 
suspension 

Figures  37  and  38  show  a  comparison  of  frequency  responses  of  both  primary 
and  secondary  suspension.  With  primary  suspension  feedback  control,  the 
results  are  the  same  as  in  Figs.  35  and  36. 

It  is  noted  that  the  active  primary  suspension  system  does  not  damp  the 
excessive  overshoots  of  secondary  suspension  in  the  transient  response  (see 
Fig.  36)  and  the  frequency  response  (see  Fig.  38).  To  achieve  the  desirable  values 
of  the  overshoots  and  the  setting  time,  a  semiactive  control  is  introduced  into  the 
secondary  suspension. 

3.6  Semiactive  Control  of  Secondary  Suspension 

Semiactive  suspension  controls  use  an  active  damper  in  parallel  with  a 
passive  spring  (see  Fig.  39).  Desired  forces  in  the  damper  are  generated  by 
modulation  via  fluid-flow  orifices  (hydraulic  or  pneumatic)  based  on  a  control 
scheme  that  uses  feedback  variables.  With  the  spring-damper  arrangement  in 
Fig.  39,  the  damper  can  produce  a  force  that  is  a  function  of  relative  velocity,  in 
our  case,  ys-yp,  of  two  masses.  What  is  required  for  significantly  improved 
control  of  transmissibility  is  a  damper  producing  a  force  proportional  to  the 
absolute  velocity,  ys,  of  the  mass.  In  physical  terms,  this  can  be  interpreted  as  the 
provision  of  a  dashpot  between  the  mass  and  an  inertial  reference,  as  shown  in 
Fig.  40.  This  arrangement,  which  has  been  described  as  a  “skyhook”  damper, 
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Fig.  37.  Frequency  response  Gp(joi)  of  primary 
suspension  using  active  feedback 
controls  in  primary  suspension 


0.1  1.0  10.0 
Frequency,  Hz 

Fig.  38.  Frequency  response  Gs(jco)  of  sec¬ 
ondary  suspension  using  active  feed¬ 
back  controls  in  primary  suspension 

offers  a  considerable  advantage  in  terms  of  transmissibility  control.  It  can  be 
constructed,  without  the  need  for  an  inertial  reference,  using  an  active  element 
under  feedback  control.  It  is  straightforward  to  show  that  this  can  be  done  by  a 
proportional  control  law  involving  the  absolute  velocity  of  the  mass.  It  was  proved 
that  semiactive  control  with  the  skyhook  configuration  can  increase  the  damping 
factor  and  that  the  resonant  peak  is  suppressed  while  high-frequency 
transmission  is  reduced  simultaneously  (Stanway  et  al.  1989). 
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Fig.  39.  Semiactive  control  for  secondary  suspensions 


Fig.  40.  Skyhook  damper  for  secondary  suspensions 
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Fig.  41.  ER  fluids  damper  based  on  shear¬ 
mode  principle 

Therefore,  the  semiactive  control  with  the  skyhook  configuration  is  desirable 
for  the  secondary  suspension  to  provide  good  ride  comfort  at  lower  costs.  This 
control  design  can  be  realized  by  a  semiactive  ER  fluid  damper  based  on  the 
shear-mode  principle  vote  Fig.  41).  In  this  configuration,  one  of  a  pair  of 
electrodes  is  fixed  to  the  piston  and  the  other  is  attached  to  the  damper  casing. 
The  relative  movement  of  the  electrodes  results  in  shear  stress  on  the  surfaces. 
The  damping  force  is  the  product  of  the  shear  stress  and  the  surface  area  of  the 
moving  electrodes.  The  shear  stress  and  hence  the  damping  force  can  be 
controlled  by  the  voltage  applied  to  the  electrodes. 

Based  on  the  principle  of  semiactive  control  of  the  secondary  suspension,  a 
feedback  control  path  kss  (active  damping)  is  added  to  the  secondary  suspension 
(see  Fig.  32).  The  root-locus  for  the  secondary  suspension  is  given  in  Fig.  42.  New 
eigenvalues  are  -3.014  ±4.020j,  with  ks  =  8.8  x  104N  s  when  using  pole  assignment 
in  the  secondary  suspension  control  design.  The  damping  ratio  of  the  secondary 
suspension  becomes  0.60. 

The  PSD  of  vehicle  acceleration  with  both  primary  and  secondary  feedback 
control  is  shown  in  Fig.  43.  Ride  comfort  is  much  improved. 


PSD  of  acceleration,  c^/H 


Fig.  42.  Root-locus  diagram  for  secondary 
suspension  feedback  control  design 


Frequency,  Hz 

Fig.  43.  PSD  of  vehicle  acceleration  using 
active  and  semiactive  feedback 
controls  in  primary  and  secondary 
suspensions  with  vehicle  speed  v  = 
100  m/s  and  guideway  roughness 
amplitude  A  =  10  6  m 
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Figures  44  and  45  show  the  comparison  of  transient  responses  of  primary 
and  secondary  suspensions  with  unit-step  input  of  guideway  disturbance. 
Response  of  the  secondary  suspension  is  improved,  and  the  setting  time  of  ts 
declines  to  0.964  s  from  2.000  s  and  overshooting  Mp  becomes  12.3%,  a  significant 
drop  from  the  original  46.3%  (see  Fig.  45  and  Table  5). 

Figures  46  and  47  show  comparisons  of  frequency  responses  of  both  primary 
and  secondary  suspension  with  both  primary  ana  secondary  feedback  control. 
These  figures  show  effects  in  reducing  vehicle  response  similar  to  those  in  Figs. 
44  and  45. 


3.7  Parametric  Analysis  of  Controi  Designs 

When  the  mechanical  control  method  is  exploited  in  maglev  suspension 
systems,  system  parameters  of  feedback  control  are  of  particular  importance  in 
depressing  oscillation  and  providing  acceptable  ride  comfort.  As  shown  in  the 
root-locus  plots  of  primary  and  secondary  suspensions  (Figs.  33  and  34),  as 
feedback  gains  increase,  system  dampings  will  decrease  and  natural  frequencies 
will  increase.  However,  as  the  feedback  gains  exceed  certain  values,  an  extra 
power  source  is  required  and  may  bring  unnecessary  cost  to  a  maglev  system. 

Although  we  do  not  extensively  discuss  the  effects  of  system  parameters  in 
control  designs,  which  we  believe  are  essential  to  any  conceptual  design  of  the 
maglev  system,  several  comparisons  on  parametric  analysis  for  control  design  of 
maglev  suspension  systems  are  performed  in  this  study. 

Figures  48-50  show  the  transient  frequency  responses  of  the  primary 
suspension  and  the  PSD  of  vehicle  acceleration,  with  various  primary  feedback 
gains  kp.  When  kp  >  1.3  x  108N,  the  oscillation  of  primary  suspension  is 
depressed.  However,  only  when  kp  >  2.7  x  108N  does  the  PSD  of  vehicle 
acceleration  reach  the  UTACV  ride  comfort  criterion. 

Figures  51-53  show  the  transient  and  frequency  responses  of  the  secondary 
suspension  and  the  PSD  of  vehicle  acceleration,  with  various  secondary  feedback 
gains  ks.  When  ks  >  6.0  x  104N  s,  setting  times  and  oscillation  depression  seem  to 
be  acceptable. 


56 


Fig.  44.  Transient  response  yp  of  primary 
suspension  with  unit-step  input  of 
guideway  perturbation  using  active 
and  semiactive  feedback  controls  in 
primary  and  secondary  suspensions 


Fig.  45.  Transient  response  ys  of  secondary 
suspension  with  unit-step  input  of 
guideway  perturbation  using  active 
and  semiactive  feedback  controls  in 
primary  and  secondary  suspensions 
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Fig.  46.  Frequency  response  Gp(jco)  of  primary 
suspension  using  active  and  semi¬ 
active  feedback  controls  in  primary 
and  secondary  suspensions 
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Fig.  47.  Frequency  response  Gs(jco)  of  second¬ 
ary  suspension  using  active  and  semi¬ 
active  feedback  controls  in  primary 
and  secondary  suspensions 
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Fig.  48.  Comparison  of  transient  response 
yp  of  primary  suspension  with  unit- 
step  input  of  guide  way  perturbation 
using  active  and  semiactive  feed¬ 
back  controls  in  primary  and 
secondary  suspensions  when  ks  = 
8.8  x  104  Ns  and  kp  =  1.3  x  108,  2.0  x 
10 8  2.7  x  108,  and  3.7  xl&N 
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Fig.  49.  Comparison  of  frequency  response  GD(jco) 
of  primary  suspension  using  active  and 
semiactive  feedback  controls  in  primary 
and  secondary  suspensions  when  ks  = 
8.8  x  104  N  s  and  kp  =  1.3  x  108,  2.0  x  108, 
2.7  x  108,  and  3.7  x  108 N 
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Frequency,  Hz 

Fig.  50.  Comparison  of  PSD  of  vehicle  acceler¬ 
ation  using  active  and  semiactive  feed¬ 
back  controls  in  primary  and  second¬ 
ary  suspensions  with  vehicle  speed  v  = 
100  mis  and  guideway  roughness 
amplitude  A  =  10  6  m  when  ks  -  8.8  x 
104  N-s  and  kp  =  1.3  x  10 8,  2.0  x  108, 
2.7  x10s,  and  3.7  x  108 N 


Fig.  51.  Comparison  of  transient  response 
ys  of  secondary  suspension  with 
unit-step  input  of  guideway  pertur¬ 
bation  using  active  and  semiactive 
feedback  controls  in  primary  and 
secondary  suspensions  when  kp  = 
2.7  x  108  N  and  ks  =  3.0  x  104,  6.0  x 
104,  8.8  x  104,  and  1.5  x  10s  N-s 
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Frequency.  Hr 

Fig.  52.  Frequency  response  Gjjcu)  of  secondary 
suspension  using  active  and  semiactive 
feedback  controls  in  primary  and  second¬ 
ary  suspensions  when  kp  =  2.7  x  108  N  and 
ks  =  3.0  x  104,  6.0  x  10 4,  8.8  x  10 4,  and  1.5  x 
105  N-s 


Frequency,  Hr 


Fig.  53.  Comparison  of  PSD  of  vehicle 
acceleration  using  active  and  semi¬ 
active  feedback  controls  in  primary 
and  secondary  suspensions  with 
vehicle  speed  v  =  100  m/s  and 
guideway  roughness  amplitude  A  = 
10  6  m  when  kp  =  2.7  x  108  N  and  = 
3.0  x  10 4,  6.0  x  10 4,  8.8  x  10 4,  and  1.5  x 
105  Ns 
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4  Closing  Remarks _ 

(1) The  dynamic  interaction  model  of  a  maglev  system  with  a  multicar, 
multiload  vehicle  traveling  along  a  flexible  guideway  was  developed  in  this  study. 
It  was  verified  that  this  model  is  desirable  for  analyses  of  vehicle/guideway 
interactions  in  maglev  systems.  The  model  can  be  incorporated  into  future 
computer  codes  for  nonlinear  dynamic  analyses  of  maglev  systems. 
Incorporation  of  this  model  into  a  computer  code  that  contains  a  six-degree-of- 
freedom  rigid-vehicle  body  has  recently  been  performed  at  Argonne  National 
Laboratory.  The  model  should  have  a  bright  future  with  many  applications  in 
commercial  maglev  systems. 

(2)  A  distributed-load  vehicle  model  is  better  than  a  concentrated-load  model 
that  may  result  in  large  amplitudes  of  both  guideway  deflections  and  vehicle 
accelerations  in  simulations.  Multicar  vehicles  have  less  car-body  acceleration 
than  does  a  single-car  vehicle  because  of  intercar  constraints.  This  indicates  that 
the  multicar  vehicle  would  provide  better  ride  comfort. 

(3)  Active  and  semiactive  feedback  control  designs  in  primary  and  secondary 
suspensions  can  be  realized  through  electro-hydraulic  systems.  The  conceptual 
designs  of  hydraulic  controllers  will  be  taken  into  account  in  our  future  work. 
The  ER  fluids  damper  may  offer  an  elegant  solution  to  the  problem  of 
implementing  semiactive  control  in  the  secondary  suspension,  which  shows 
great  potential  for  reducing  costs  and  improving  ride  comfort. 

(4)  Maglev  may  become  a  major  transportation  mode  in  the  21st  century. 
Because  the  cost  for  a  commercial  maglev  system  is  still  very  high,  it  is  wise  to 
consider  dynamic  control  systems  before  completing  the  guideway  design  so  that 
overall  system  cost  can  be  reduced. 
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